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A canonical partition function for the two-component excluded volume model is derived, leading 
to two different van der Waals approximations. The one is known as the Lorentz-Berthelot mixture 
and the other has been proposed recently. Both models are analysed in the canonical and grand 
canonical ensemble. In comparison with the one-component van der Waals excluded volume model 
the suppression of particle densities is reduced in these two-component formulations, but in two 
essentially different ways. Presently used multi-component models have no such reduction. They 
are shown to be not correct when used for components with different hard-core radii. 
For high temperatures the excluded volume interaction is refined by accounting for the Lorentz 
contraction of the spherical excluded volumes, which leads to a distinct enhancement of lighter 
particles. The resulting effects on pion yield ratios are studied for AGS and SPS data. 
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I. INTRODUCTION 

Thermal models are commonly used to interprete ex- 
perimental data from hadron collisions, see for instance 
[1-5] . In the van der Waals excluded volume (VdW) 
model the short range repulsion between particles is rep- 
resented by hard-core potentials, i. e. the finite size of the 
particles is taken into account. As a consequence par- 
ticle yields are essentially reduced in comparison with 
ideal gas results, whereas yield ratios remain almost un- 
changed, if the same hard-core radius is used for all par- 
ticle species. 

As particle species with smaller hard-core radii are 
closer to the ideal case, their particle densities are sup- 
pressed less. Consequently, their yield ratios to particle 
species with larger hard-core radii are enhanced. This 
fact has been used in recent efforts [6] to explain the 
experimentally observed pion abundance for AGS and 
SPS data [7] by introducing a smaller hard-core radius 
for pions R w than for all other hadrons R . However, 
the resulting values are quite large, R = 0.8 fm and 
= 0.62 fm. In Ref. [8] a reasonable fit of SPS data 
has been obtained only for a distinctly smaller pair of 
hard-core radii. 

The excluded volume models used in [6,8] , however, 
are not correct in the case of different hard-core radii. 
As will be shown in Sect. II , these models correspond to 
a system where the components are separated from each 
other by a mobile wall and hence cannot mix. 
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A more realistic approach requires a two-component 
partition function including a term for the repulsion 
between particles of different hard-core radii. In the 
case of two components, however, the VdW approxima- 
tion is not uniquely defined. The simplest possibility 
yields the Lorentz-Bcrthclot mixture, which was origi- 
nally postulated by van der Waals for binary mixtures, 
see Refs. [9-12] . Another VdW approximation was re- 
cently proposed in Ref. [13] . These two formulations 
contain a suppression of particle densities similar to the 
one-component van der Waals gas, which is reduced to 
different extend for each formulation. In the present work 
we will study and apply both of these formulations. 

There is yet another cause for a reduced excluded vol- 
ume suppression. Particles are considered to be rigid 
spheres in the VdW model. At high energies as achieved 
in nuclear collisions, however, relativistic effects cannot 
be neglected [14] . Within the logic of the VdW model 
it is necessary to take into account the Lorentz contrac- 
tion of the spheres. We will use an approach developed 
in Ref. [15] providing approximative formulae for rela- 
tivistic excluded volumes: naturally, they decrease with 
rising temperature, and the effect is stronger for lighter 
particles. At high temperatures, consequently, it is not 
possible to use a one-component VdW description (i. e. a 
common excluded volume for all particle species) for a 
system of species with various masses. Since different 
masses cause different reductions of the excluded vol- 
umes at a given temperature, a multi- component VdW 
description is required. 

To illustrate the influence of different excluded volumes 
we will restrict ourselves in this work to the simplest 
'multi-component' case, the two-component case. The 



1 



crucial extension from the one- to the two-component 
case is to include the repulsion between particles of two 
different hard-core radii. As it will be illustrated, a gen- 
eralisation to the multi-component case is straightfor- 
ward and will yield no essential differences [16] . 

In the next section a derivation of the one-component 
canonical partition function with (constant) excluded 
volumes is presented. The generalisation to the two- 
component case is made and two possible VdW approxi- 
mations are analysed: the Lorcntz-Bcrthclot mixture [12] 
and the recently proposed approximation of Ref. [13] . 
The corresponding formulae for the grand canonical en- 
semble are derived and discussed in Sect. Ill . Relativis- 
tic excluded volumes are introduced in Sect. IV and the 
corresponding equations of state are presented for both 
models. In Sect. V a fit of particle yield ratios [6] is 
re-evaluated for both approximations, with constant and 
with relativistic excluded volumes. The conclusions are 
given in Sect. VI . 

The appendixes A and C give a detailed analysis and 
comparison of the two approximations in the canonical 
and grand canonical ensemble, respectively. In App. B 
a general proof of the thermodynamical stability of the 
non-linear approximation is given. 



II. CANONICAL TREATMENT 

First we will derive the canonical partition function 
(CPF) for the one-component VdW gas by estimating 
the excluded volumes of particle clusters. Then this pro- 
cedure will be generalised to the two-component case. 

For simplicity Boltzmann statistics are used through- 
out this work. The deviations from quantum statistics 
are negligible as long as the density over temperature 
ratio is small. This is the case for the hadron gas at tem- 
peratures and densities typical for heavy ion collisions, 
see e. g. Ref. [6] . 

Note that in this work we will use the term VdW for 
the van der Waals excluded volume model, not for the 
general van der Waals model which includes attraction. 



A. The Van der Waals Excluded Volume Model 

Let us consider N identical particles with temperature 
T kept in a sufficiently large volume V , so that finite 
volume effects can be neglected. The partition function 
of this system (U = c = fcs = 1) reads 

6 N r 

Z(T,V,N) = ^ y d 3 *i ••• d 3 x N cxp[-^] . (1) 

Here, <j> = <fi(T; to, g) denotes the momentum integral of 
the one particle partition 



4>(T; m, g) = ^ / dfc fc 2 exp 



E(k) 
T 



where E{k) = \/k 2 + to 2 is the relativistic energy and 
g = (25 + 1)(2J + 1) counts the spin and isospin de- 
generacy. For a hard-core potential Un of N spherical 
particles with radii R the potential term in Eq. (1) reads 



exp[-^] - II 0{\x l3 \-2R) , 



(3) 



i<j<N 



where xy denotes the relative position vector connecting 
the centers of the i-th and j-th particle. Hence one can 
write 

/ d 3 xi ••• d 3 x N exp[-^l 
Jv N 

= [ d 3 Xl ••• d 3 x N J] 0(l*«|-2fl) 

JV N l<i<j<N 

= [ d 3 Xl [ d 3 x 2 8(\x 12 \ - 2R) x 
Jv Jv 

x ••• / d 3 x N n 0(\x itN \-2R) 

JV 1<»<JV-1 

{Si} {xi...x N _i} 

= J d 3 Xl J d 3 X 2 ••• J d 



XN ■ 



(4) 



{Si...Xj} 

Here, J d 3 Xj+i denotes the available volume for 
Xj+i, which is the center of the particle with num- 
ber j + 1 , if the j other particles are configurated as 
{xi . . . Xj} . We will show now that this volume is es- 

{x!...Xj} 

timated by J d 3 x J+ i > (V — 2b j) , where 2b = 



(2) 



^-(2R) 3 is the excluded volume of an isolated particle 
seen by a second one. Then, 2b j estimates the total vol- 
ume which is excluded by all particle clusters occuring in 
the configuration {x*i . . . Xj} . 

It is sufficient to prove that the excluded volume of 
a cluster of k particles is less than the excluded volume 
of k isolated particles. A group of k particles forms a 
fc-cluster, if for any of these particles there is a neigh- 
bouring particle of this group at a distance less than AR . 
The exact excluded volume of a fc-cluster, V(u\ , obvi- 
ously depends on the configuration of the k particles. If 
one considers two isolated particles, i. e. two 1-clusters, 
and reduces their distance below 4i? , their excluded vol- 
umes will overlap. They form now a 2-clustcr with the 
excluded volume «( 2 ) = 4b— lv ov , where v ov denotes the 
overlap volume. 

Evidently, one can construct any fc-cluster by attaching 
additional particles and calculate its excluded volume by 
subtracting each occuring overlap volume from 2b k . It 
follows that V(u\ < 2b k is valid for any fc-cluster, and 
this inequality leads to the above estimate. Obviously, 
its accuracy improves with the dilutcness of the gas. 

Using these considerations one can approximate the 
r. h. s. of Eq. (4) : starting with j + 1 = N one gradually 
{xi...Sj} 

replaces all integrals J d 3 xj + i by (V — 2b j) . One 
has to proceed from the right to the left, because only the 
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respective rightmost of these integrals can be estimated 
in the described way. Hence one finds 



lN N-l 



Z(T,V,N)>^- J] (V-2bj). 



(5) 



In this treatment the VdW approximation consists of 
two assumptions concerning Eq. (5) . Firstly, the product 
can be approximated by 



The pressure formula (11) corresponds to the Boltz- 
manii approximation of the commonly used two- 
component VdW models of Refs. [6,8] as will be shown 
in Sect. V. It is evident that p sp (11) does not describe 
the general two-component situation without a separat- 
ing wall. Therefore, it is necessary to find a more realistic 
model, i. e. an approximation from a two-component par- 
tition function. This will be done in the following. 



JV-l 



n (i-ti) 

3=0 



exp 



N-l 2b ■ 



EiV — 
3=0 



V 



(6) 



= exp [_£(JV-l)JV]£*(l-£iV) 



N 



where exp [—a;] = (1 — x) is used for dilute systems, i. c. 
for low densities 2bN/V <C 1 . The second assumption is 
to take the equality instead of the inequality in Eq. (5) . 
Then the CPF takes the VdW form, 



Z vdw (T,V,N) 



N 



N\ 



(V - bN) 



N 



(7) 



As usual, the VdW CPF is obtained as an approximation 
for dilute systems, but when used for high densities it 
should be considered as an extrapolation. 

Finally, one obtains the well-known VdW pressure for- 
mula from the thermodynamical identity p(T, V, N) = 
Tdha[Z(T,V,N)]/dV, 



V vdw -(T,V,N) 



TN 
V -bN 



(8) 



using the logarithm of the Stirling formula. 

Now let us briefly investigate a system of volume V 
containing two components with different hard-core radii 
R\ and R 2 which are separated by a wall and occupy the 
volume fractions xV and (1 — x)V , respectively. Accord- 
ing to Eq. (8) their pressures read 



PvdwCT.zV.JVi) = 
Pvdw(T,(l-x)V,N 2 ) = 



T N\ 



xV-Nrfu ' 
(1 - x)V - N 2 b 2 



(9) 
(10) 



where the particle numbers Ni , N 2 and the excluded vol- 
umes 6n = i|2L _r3 j 5 22 = 1|zl _r| correS p 0n d to the 
components 1 and 2 , respectively. 

If the separating wall is mobile, the pressures (9) and 
(10) must be equal. In this case the fraction x can be 
eliminated and one obtains the common pressure of the 
whole system 



Pvdw(T, xV, Ni) = pv dW (T, (1 - x)V, N 2 ) 

N 2 ) 



= P sp (T,v,m,N 2 )= T{Nl 



V - Nibu - N 2 b 



22 



(11) 



Since the components are separated in this model system 
it will be referred to as the separated model [16] . 



B. Generalisation to the Two-component Case 

Recall the simple estimate (4-7) , which gives a phys- 
ically transparent derivation of the one-component CPF 
in the VdW approximation. Let us use it now for a 
two- component gas of spherical particles with radii R\ 
and R 2 , respectively It is important to mention that 
each component may consist of several particle species 
as long as these species have one common hard-core ra- 
dius, i. e. the number of necessary VdW components is 
determined by the number of different excluded volume 
terms b qq . In the case of two radii the potential term (3) 
becomes 



exp 



Un 1 + n 2 
T 



= n 0{\x ij \-2R 1 )x (12) 

i<j<Ni 

x J] d{\x u \-2R 2 ) x 

k<e<N 2 

x n 0(\x mn \ - (Ri + R 2 )) . 

m<Nx 
n<N 2 



The integration is carried out in the way described above; 
e. g. firstly over the coordinates of the particles of the sec- 
ond component, then over those of the first component. 
For the estimation of the excluded volume of a /c-cluster 
now two different particle sizes have to be considered. 
One obtains 



Z(T,V,N lt N 2 ) 

> 



n (^-2611*)^ x 

i=0 



N^. N 2 \ 
fjv 2 -i 

x { n (^-2612^1-2622^ 
3=0 



^ <H <P 2 yWi+A^ 

JVi! N 2 l 



x exp 



jV 1 2 bn+2Af 1 jV 2 b 12 +jV 2 2 b 22 
V 



(13) 



where it is (f> q = (j)(T; m q , g q ) , and 2b pq = ^ (R p + R q ) 3 
denotes the excluded volume of a particle of the compo- 
nent p seen by a particle of the component q (p, q = 
1,2 hereafter). Approximating the above exponent by 
exp [—a;] = (1 — x) yields the simplest possibility of a 
VdW approximation for the two-component CPF , 

Z^ dw (T,V,N u N 2 ) 



3 



K N1 ^ 



N\\ N 2 l 



(14) 



x V 



Nfbu + 2N 1 N 2 b 12 + N 2 2 b 22 

m + N 2 



Nil N 2 l 

( NN 2 
x I V- Nibn - N 2 b 22 + - t D 



(15) 



N x +N 2 



Ni+N 2 

where the non-negative coefficient D is given by 



D = bu + b 22 - 2 b 



12 



(16) 



This approximation will be called the non-linear approx- 
imation as the volume correction in (15) contains non- 
linear terms in N\ , N 2 . The corresponding pressure fol- 
lows from the thcrmodynamical identity, 



p n \T,V,Ni,N 2 ) 



Pi'+P 2 



(17) 



T{Ni+N 2 ) 



V - Nibu - N 2 b 22 + 



Ni N 2 
N ± +N 2 



D 



This canonical formula corresponds to the Lorentz- 
Bcrthclot mixture (without attraction terms) known 
from the theory of fluids [12] . ft was postulated by van 
der Waals [9] and studied as well by Lorentz [10] and 
Berthelot [11] . 

The crucial step from the one- to the two-component 
gas is to include b pq terms (p ^ q) additionally to the 
b qq = b\ R=Rq terms. For the multi-component gas no 
further essential extension is necessary. Consequently, 
the generalisation of the above procedure to the multi- 
component case, i. e. an arbitrary number of different 
hard-core radii, is straightforward [16] . 

In Ref. [13] a more involved approximation has been 
suggested for the two-component VdW gas. This follows 
from splitting the exponent in the CPF (13) by intro- 
ducing two generalised excluded volume terms bi 2 and 
b 2 i (instead of a single and symmetric term 2 612) for the 
mixed case, 



Z(T,V,N U N 2 ) 

1 Ni ,N 2 

^ fi_ <h_ 

~ Nil N 2 l 



V 



N x +N 2 



(18) 



x exp 



N^b 11 +N 1 N 2 (b 12 +b 21 )+N^b 22 
V 



which leads to an alternative two-component VdW CPF , 



^ww {T, V, Ni,N 2 ) 

= 

~ Ni 



i-j- (V - Nibu - A 2 6 2 i) 1 x 
(V - N 2 b 22 - JV1612) 



(19) 



\ N 2 



N 2 l 



Since the particle numbers Ni , N 2 appear solely linearly 
in the volume corrections, these formulae will be referred 



to as the linear approximation. In this approximation 
one obtains [13] for the pressure 



p Vln (T, V, Ni,N 2 ) 
TNi 



Ain 



-v 2 



(20) 



■ + 



TNo 



V - Nibu - N 2 b 2 i V - N 2 b 22 - Nibi 2 



The choice of the generalised excluded volume terms 
bpq is not unique in the sense that all choices which satisfy 
the basic constraint 612 + & 2 i = 2&12 are consistent with 
the second order virial expansion [13] . Therefore, ad- 
ditional conditions are necessary to fix these generalised 
excluded volumes. In Ref. [13] they were chosen as 



b 12 = bu 



2 bi 



bu + b 



hi = b 2 



2 bi 



'22 



hi + 622 



(21) 



For this choice, the linear approximation reproduces a 
traditional VdW gas behaviour, i. e. one- component-like, 
in the two limits R 2 — Ri and R 2 — as readily checked. 
The factor 2&i 2 /(&n + b 22 ) = 1 - D/(hi + b 22 ) is al- 
ways smaller than unity for Ri ^ R 2 , consequently, the 
b pq terms are smaller than the corresponding terms b pp . 
Note that there are many possible choices for h 2 and 
621 , e. g. additionally dependent on the particle numbers 
Ni and N 2 , whereas the non-linear approximation (14) 
contains no such additional parameters. 

The formulae of the linear approximation are gener- 
ally valid for any choice of h 2 and 621 satisfying the con- 
straint 612 + hi = 2b i2 . In the following, however, we 
will restrict our study to the special choice given in the 
Eqs. (21) . The canonical (and grand canonical) formulae 
for the multi-component case are given in Ref. [13] . 



C. Comparison of both 
Two-component VdW Approximations 

As the VdW approximation is a low density approxi- 
mation it is evident that the linear and non-linear formu- 
lae are equivalent for such densities. Deviations, however, 
occur at high densities, where any VdW approximation 
generally becomes inadequate. 

The differences between both approximations result 
from the fact that the linear pressure (20) has two poles, 
Vi — V and v 2 n — V , whereas the non-linear pressure 
(17) has solely one pole, v nl — V . For constant vol- 
ume V these poles define limiting densities, e. g. hi = 
max(7Vi/y) as functions of n 2 = N 2 /V , 



. lin 



(NuN 2 ) = V hi{n 2 ) = h^{n 2 ) (22) 

or v n \Ni,N 2 ) = V ~» hi(n 2 ) = n^) , (23) 

which represent the domains of the two pressure formulae 
in the 712-ni-plane. The explicit fomulae are discussed 
in App. A . 

In Fig. 1 (a) an example of these limiting densities is 
shown for R 2 /Ri = 0.4 . It is clearly seen that the non- 
linear domain (below the solid line) is larger than the 
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linear domain (below both dashed lines), which is gen- 
erally the case for R 2 ^ R\ . Especially for R 2 <C R\ 
the non-linear domain is distinctly larger for high den- 
sities of the large component, n\b\\ > n 2 b 22 , whereas 
both domains are similar for high densities of the small 
component, n 2 b 2 2 > ^l&n ■ 

The linear approximation is constructed in traditional 
VdW spirit; the densities n q ln achieved in this approx- 
imation are below the maximum density of the corre- 
sponding one- component VdW gas max(n° c ) = l/b qq , 
which is defined by the pole of p° c = Pvdw(T, V, N q ; b qq ) 
from Eq. (8) . 

In the non-linear approximation, however, the possi- 
ble densities of the larger particles n" 1 can exceed 1/fen 
due to the occurence of negative partial derivatives of 
the pressure, dp nl /dN 2 < 0. In this context it is nec- 
essary to state that this behaviour does not lead to a 
thermodynamical instability of the non-linear approxi- 
mation as proven in App. B . The linear approximation 
shows no such behaviour, it is always dp hn /dNi > and 
dp 1[n /dN 2 > 0. 

The condition dp nl /dN 2 = defines the boundary 
n^ l ' hd (n 2 ) of the region of negative partial derivatives 
of the non-linear pressure. In Fig. 1 (a) this boundary is 
shown by the dotted line for R 2 /R\ = 0.4; the values of 
dp nl /dN 2 are negative above this line. 

Densities larger than n" 1 = 1/bn can only occur, if R 2 
is smaller than a critical radius, 

R2 < i?2,critORl) = (#4- l)Ri « -Ri/1-7 . (24) 

Then, the boundary n^' hd (n 2 ) starts inside the non- 
linear domain, see App. A for details. 

The reason for this behaviour is the ratio of the 
amounts of small and large particles. There are much 
more small than large particles in the system for densi- 
ties ni,7i2 along the boundary ' (n 2 ) at high den- 
sities ri\ : here, the fewer large particles are surrounded 
by many small particles. Therefore, the excluded volume 
interaction of the large particles in the non-linear pres- 
sure (17) is governed not by the simple term &n but by 
the mixed term 612 , which is distinctly smaller than b\\ 
for R 2 -C i?i . The maximum density achieved in the 
non-linear approximation max(nf') = 4/&n is obtained 
for i? 2 — > and N 2 » N\ , i. e. these formulae go far be- 
yond the traditional VdW results in the corresponding 
situation. 

An example of pressure profiles for p\ m , p 2 n and p nl 
for ni&H = 0.9 is shown in Fig. 1 (b) , where it is 
R 2 /Ri — 0.4 as in Fig. 1 (a) . The non-linear pressure 
(solid line) firstly decreases as the densities ni,n 2 cor- 
respond to the region of negative partial derivatives, see 
Fig. 1 (a) . The partial pressures of the linear approxima- 
tion are shown by dashed lines. The non- linear domain 
is seen to be larger since it is one of the linear partial 
pressures which diverges first for increasing n 2 . 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
n 2 ■ b 22 

(b) 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 
n 2 • b 22 



FIG. 1. (a) Domains of the linear and non-linear approx- 
imation for R2/R1 = 0.4: limiting densities hi (isobars for 
p q (ni, n'z) = 00) and the lower boundary h™ l ' hd to the region 
of negative partial derivatives of the non-linear pressure. The 
dashed lines correspond to the two poles of the linear pressure, 
and the solid line corresponds to the pole of the non-linear 
pressure. For given U2 the possible densities nl ln are below 
both dashed lines, whereas the possible densities nf are be- 
low the solid line. Negative derivatives dp nl /dN2 < occur 
only above the dotted line. 

(b) Pressure profiles in dimensionless units for R2/R1 = 0.4 as 
in (a) at fixed ni&n = 0.9 . The dashed lines show the partial 
pressures of the linear approximation p\ ln and p^ 111 , while the 
solid line shows the total pressure of the non-linear approxi- 
mation p ul with initial decrease due to negative dp nl /dN2 . 
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We conclude that the linear and non-linear approxi- 
mation show a drastically different behaviour for high 
values of the large component's density ri\ . In the lin- 
ear approximation (20) the possible density values are 
below l/6n and l/b 22 , respectively, and the derivatives 
dp hn /dN q are always positive. Whereas in the non- 
linear approximation (17) higher densities m > I/&11 
are possible due to the occurence of negative derivatives 
dp nl /dN 2 < 0. This may be considered as pathologi- 
cal - or used as an advantage to describe special situa- 
tions, e.g. densities 1/bn < n\ < n" 1 for R 2 -C Ri (see 
App. A). 

However, the use of any Vd W approximation is in prin- 
ciple problematic for densities near l/b qq . For low den- 
sities the non-linear and linear approximation are prac- 
tically equivalent, and the non- linear approximation is 
preferable since the formulae are essentially simpler. 



where A q = A(T, fi q ; m q , g q ) = exp[^ g /T] 4> q . 

The evaluation of both maximum conditions for the 
grand canonical pressure 







d 



In 



A JVi , M 2 
^1 /i 2 

aTTT "7V2T 



(28) 



N1+N2 



yields a system of two coupled transcendental equations, 



= A\ exp 



-(€i nl + C 2 nl )fon + 



«r+«2' 



D 



e 2 nl (T ;Ai i, M2 ) 



A 2 exp -{S? + S?)b22 



D 



(29) 
(30) 



III. GRAND CANONICAL TREATMENT 

Let us now turn to the grand canonical ensemble 
(GCE) . The grand canonical partition function is built 
using the CPF , 



OO OO 

= E E exp 

AT 1= N 2 =0 



(25) 



T 



Z(T,V,N U N 2 ) 



whereas the chemical potentials [i\ and [i 2 correspond to 
the components 1 and 2 , respectively. For the VdW CPF 
(15) or (19) there are limiting particle numbers N\(N 2 ) or 
N 2 (Ni) , where each CPF becomes zero. For this reason 
the above sum contains only a finite number of terms. 
Then it can be shown that in the thermodynamical limit 
(i. e. in the limit V —* 00 for N q /V = const.) the grand 
canonical pressure p(T, fj,i, fi 2 ) = T \n[Z{T 1 V, Hi, fj, 2 )]/V 
depends only on the maximum term of the double sum 
(25) , where Ni = M and N 2 = Mi . One obtains 



P(T, (J-i, IJ-2) 



(26) 



T 

= lim — In 

y^oo V 



exp 



T 



z{t,vMiM2) 



wheras N\ and Af 2 are the average particle numbers. 



A. The Two VdW Approximations 

For the non-linear VdW approximation (15) the last 
expression takes the form 



p n \T^ llt i 2 ) 
T 

= lim — In 

V^OD V 



(27) 



U7T "a^T 



x (v -N^-N^ + ^^d) 



M+AT 2 



where and £ 2 arc defined as 



£1 
£2" 



ni 



M 



V-M&n -N 2 b 22 
M 2 



M X N 2 n 



V-Mbn-Mb^ + ^^D 



(31) 
(32) 



In the thermodynamical limit the average particle num- 
bers Afi and N 2 are proportional to V as N q — n™ 1 V . 
Then the volume V disappears in the definitions of £f l 
and £ 2 nl given by Eqs. (31, 32) , and they can be solved 
for either the density n" 1 or n 2 l , 



»i {T, fJ-i, fJ-2) = 



n 2 \T,m,Li 2 ) = 



£1 



nl 



tnl 
S2 



1 + £i%i + U il b22 - fr^r D 



(33) 
(34) 



The £ g nl = £g nl (T, /ii, fj, 2 ) are the solutions of the coupled 
Eqs. (29) and (30) , respectively. 

Hence, the pressure (27) can be rewritten in terms of 
£i nl (29) and £ 2 nl (30) , 



p nl (T, Ml , M2 )=T (£i nl + £ 2 nl ) 



(35) 



supposed that Eqs. (33, 34) are taken into account. If 
the definitions (31) and (32) are used for and £ 2 nl , 
the pressure formula (35) coincides with the canonical 
expression (17) for Ni — N\ and A 2 = J\f 2 . 

Since the formulation is thermodynamically self- 
consistent the identity n q = dp(T, fi\, fj, 2 )/dfj, q leads to 
Eqs. (33, 34) as well. 

The grand canonical formulae of the linear approxima- 
tion [13] are obtained exactly as presented for the non- 
linear case in Eqs. (27-35) . In the linear case the two 
coupled transcendental equations are 



G 



e 1 lin (T, Ml , M2 ) = A 1 exp 



£ 2 lin (T, m , M2 ) = A 2 cxp [-£ 2 lin 6 22 - 621 
with the definitions 



lin 



M 



V - M&11 - A/" 2 6 2 i 



■lin 



7V 2 



V - Af 2 b 22 -Wife 



(36) 
(37) 

(38) 
(39) 
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The expressions for the particle densities are found by 
solving Eqs. (38, 39) for either n{ in or , 



nl in (T, Ml , M2 ) 



(40) 



£i lin (l+£ 2 lin [622-621]) 



1 + £i lin 6n + e 2 lin 6 22 + £i in £ 2 lin [6n6 22 - 612621] ' 
nf(T, Ml , M2 ) (41) 

= e 2 lin (i+ei lin [611-612]) 

1 + £i Iin 6n + ^""622 + ei lin e 2 lin [611622 - 612621] ' 

For the linear approximation the pressure (26) can be 
rewritten in terms of £f ln (36) and £ 2 lln (37) , 



p lin cz>i, M2 ) = T (er + eJ 



(42) 



if Eqs. (40, 41) are taken into account, like in the non- 
linear case. 



B. Comparison of both Approximations 

Let us briefly return to the usual VdW case, the 
one-component case. The corresponding transcendental 
equation is obtained from either Eqs. (29, 30) or (36, 37) 
by setting i?i = R 2 = R and A x = A 2 = A , 



r c (T,M) = AexpK oc 6] 



(43) 



whereas b = bu = 622 • The transcendental factor 
cxp[— £ oc b] has the form of a suppression term, and the 
solution £ oc = p oc /T of this transcendental equation ev- 
idently decreases with increasing 6 for constant T and 
/i . Then in turn, the corresponding particle density 
n oc _ £° C /(1 _|_ £° c 5) i s suppressed in comparison with 
the ideal gas due to the lower £ oc and the additional 
denominator. Thus, a suppressive transcendental factor 
corresponds to a suppression of particle densities. 

Now it can be seen from Eqs. (29, 30) and (36, 37) that 
the transcendental factors of both two-component ap- 
proximations contain as well this usual one-component- 
or VdW-like suppressive part exp[— (p/T) b qq ] . But since 
it is D > and b pq < b pp , respectively, there is further- 
more an attractive part in each corresponding transcen- 
dental factor. 

In the non-linear approximation the attractive part can 
even dominate the suppressive part for the smaller com- 
ponent, e. g. in Eq. (30) for R 2 < R\ . Then the larger 



component can reach densities n" 1 higher than l/6n, 
analogous to the CE . A detailed discussion is given in 
App. C . 

High densities in the canonical treatment correspond 
to large values of the chemical potentials in the grand 
canonical treatment. In the limit 

m/T -> 00 (T, /i 2 = const.) or -> 00 (44) 

the solution of Eq. (30) , £ 2 nl , can be enhanced for in- 
creasing instead of being suppressed, if R 2 is suffi- 
ciently small. This may be called the non-linear enhance- 
ment. The behaviour of the non-linear approximation in 
the limit (44) depends only on the ratio of the two radii 
i?i / R 2 and is characterised by the coefficient 



«2 = V D / b 22 - 1 • 



(45) 



A negative a 2 relates to a suppressive transcendental fac- 
tor in Eq. (30) . For equal radii R 2 = R\ it is a 2 = — 1 , 
and the suppression is evidently not reduced but VdW- 
like. For —1 < a 2 < this suppression is reduced, the 
most strongly for a 2 w . 

In the case a 2 = the suppression for £ 2 nl (30) disap- 
pears in the limit (44) , on has £ 2 nl — > A 2 = const. This 
case provides the critical radius R 2t cr ;t (24) . 

For a 2 > or R 2 < R 2 , cr it the non-linear enhancement 
of £ 2 nl occurs for increasing ; it is the stronger the 
larger a 2 is. Then nf 1 (33) can exceed max(n° c ) = l/6n , 
whereas n 2 l (34) does not vanish (see App. C for the 
explicit fomulae). The density max(n" 1 ) = 4/6n is 
achieved for a 2 — > 00 or R 2 — > . 

The suppression in the transcendental factor of (29) 
is generally reduced for R 2 < R\ , the more strongly the 
smaller R 2 is, but there is no enhancement possible in 
the limit (44) . 



IV. RELATIVISTIC EXCLUDED VOLUMES 

In this section we will investigate the influence of rc- 
lativistic effects on the excluded volumes of fast moving 
particles by accounting for their ellipsoidal shape due to 
Lorentz contraction. In Ref. [15] a quite simple, ultra- 
relativistic approach has been made to estimate these 
effects: instead of ellipsoids two cylinders with the corre- 
sponding radii have been used to calculate approximately 
the excluded volume term b pq for the two-component 
mixture. The resulting relativistic excluded volumes de- 
pend on the temperature and contain the radii and the 
masses as parameters. The simple, non-mixed term reads 
[14] 



/ 3771" a„ 7T 2 , 

b gg (T) = a qq ( — f + -)Rj 



(46) 



where a q = a(T;m q ,g q ) denotes the ideal gas scalar den- 
sity, 
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a(T;m,g) 



9 
2tt 2 



dfc k 2 



m 

W) 



exp 



E(k) 
T 



(47) 



The expression for the mixed case can be derived sim- 
ilarly from [15] , 

&i 2 (T) = a 12 {(^/ 1 + ^f)i? 1 3 (48) 

whereas the abbreviations /i and f 2 are dimensionless 
functions of both radii, 

A = f(2 + ^ + S*). /2 = f(2 + ^ + |f) . 



The normalisation factors 

16 



an = «22 



an 



IT 

3 



2 



7T 2 fi 2 



(49) 
(50) 



are introduced to normalise the ultra-relativistic approx- 
imations (46, 48) for T = to the corresponding non- 
relativistic results. For the hadron gas, however, these 
Boltzmann statistical formulae will only be used at high 
temperatures, where effects of quantum statistics are neg- 
ligible. 

Note that it is not appropriate to consider tempera- 
ture dependent hard-core radii R P {T) or R q (T) since the 
b pq (T) terms give the Lorente-contracted excluded vol- 
umes and are involved functions of T 7 m p ,m q7 R p and 
R q . However, for a given value of b pq (T) the necessary 
hard-core radii R p and R q will obviously depend on the 
temperature. 

It is evident that the formulae (46, 48) suffice already 
for the multi-component case, because even a multi- 
component VdW formulation contains only b pq terms. 

For both approximations the expressions for the pres- 
sure (35) or (42) and corresponding particle densities (33, 
34) or (40, 41) remain unchanged. However, due to the 
temperature dependence of the relativistic excluded vol- 
umes the entropy density is modified 

d 

= S nre ] + S re l (0T&11 , d T b 2 2 , d T h 2 ) ■ ( 51 ) 

The additional term s ro i depends on temperature deriva- 
tives of the relativistic excluded volumes, drb pq = 
db pq /dT , which represent their thermal compressibility. 

Furthermore, the term s rc i generates additional terms 
for the energy density, according to e = Ts — p + H\n\ + 
\i 2 n 2 ■ In the non-linear approximation one obtains 



(T,fn,fj, 2 ) 

= nf 5- + < £ - K 1 + n 2 nl ) T 2 x 

01 4>2 

x (V &rb n + £ 2 nl d T b 22 - d T D) 



(52) 



and the linear approximation yields 
e lin (T, m , M2 ) = nl in ^+n^ ' 2 



lin rji 2 ( /- lin o i , /■ lin 



(53) 



-n[ ln T 2 ^rd T bu+^ 2 hn d T b 12 
-n' 2 in T 2 (ttd T b 22 +^d T b 21 ) , 

whereas e q = e(T;m q7 g q ) denotes the ideal gas energy 
density 



e(T;m,g) = -^-^ j dkk 2 E(k) exp 



E(k) 
T 



(54) 



The additional terms in the entropy density (51) and in 
the energy density (52) or (53) which contain tempera- 
ture derivatives do evidently not occur in the case of the 
usual non-relativistic, i. e. constant excluded volumes. 

Let us now study the hadronic equation of state gen- 
erated by each of the two-component VdW approxima- 
tions and their modifications due to relativistic excluded 
volumes. When used to describe hadronic particles, the 
hard-core radii R q should be considered as parameters 
rather than particle radii. We identify the first compo- 
nent as nucleons (mi = m n = 939 MeV, yUi = fi n = 
and gi = g n = 4 for symmetric nuclear matter) and the 
second as pions (m 2 = = 138 MeV, fi 2 = /i^ = 
and g 2 = g v = 3) . Quantum statistical effects other 
than the degeneracy factors g q are neglected. To repro- 
duce experimental data, however, it would be necessary 
to consider all hadrons and hadronic resonances as well 
as the contributions from hadronic decays into daughter 
hadrons. 

For some examples the temperature dependence of 
the relativistic excluded volumes is shown in Fig. 2 (a) , 
given in units of the corresponding non-relativistic terms, 
b pq — b pq (0) . The solid line and the short dashes show 
the basic excluded volumes &n(T) and b 22 (T) , respec- 
tively. In these relative units the decreases of 6n(T) and 
b 22 {T) depend only on the corresponding masses. It is ap- 
parent that the pion excluded volume b 22 (T) is affected 
much stronger than the excluded volume of the nucleons, 
6n(T) . The dotted line shows the mixed volume term 
&12 (T) , and the long dashes show the compound volume 
tcrmL>(T) = bn{T) + b 22 {T)-2b 12 {T) . These two terms 
depend obviously on both masses and both radii. 

The curves for the generalised excluded volume terms 
of the linear approximation b\ 2 {T) and b 2 \{T) behave 
similarly to b\ 2 (T) . 

Introducing the relativistic excluded volumes b pq {T) , 
however, has two effects. First, the maximum densities 
become larger since it is generally l/b qq (T) > l/b qq as 
seen in Fig. 2 (a) . Furthermore, the balance between 
the lighter and the heavier species is changed because 
the lighter species is affected more than the heavier at 
the same temperature: For the above parameters it is 
b 22 (T)/b 22 <bn(T)/b u . 
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FIG. 2. Temperature dependence of the relativistic ex- 
cluded volume terms for mi = m n , m 2 = ra^,R\ = 0.6 fm . 

(a) Relative values for R 2 = 0.3 f m : 6ii(T)/6n , b 12 {T)/b 12 , 
b 22 (T)/b 22 and D(T)/D (solid line, dotted line, short and 
long dashes, respectively). The relativistic excluded volume 
of light species (b 22 (T)/b22) is affected more strongly by tem- 
perature. 

(b) The characteristic coefficient of the non-linear approxi- 
mation a 2 (T) = {yfD(T)/b 22 (T) - 1) for various R 2 = 0.25 , 
0.3 and 0.4 fm . The non-linear enhancement (a 2 (T) > 0) be- 
comes stronger due to the decrease of the relativistic excluded 
volumes with increasing temperature. 



In the non-linear approximation this balance is char- 
acterised by the coefficient a 2 defined by Eq. (45) . 
In Fig. 2(b) the temperature dependence of a 2 (T) = 
{yf 'D{T) /b 22 {T) — 1) is shown for three different values 
of R 2 • The relativistic coefficient a 2 (T) increases with 
T , i. e. the non-linear enhancement becomes stronger for 
higher temperatures. For some values of R 2 , e. g. R 2 = 

0. 4 fm, a primary suppression a 2 (0) = a 2 < 0, turns 
into an enhancement a 2 (T) > , when the temperature 
is sufficiently high. For temperature dependent excluded 
volumes i?2,crit looses its meaning; here, only a 2 (T) > 
is the valid condition for the occurrence of the non-linear 
enhancement or densities n" 1 > l/6n(T) . 

The linear coefficient, a 2 (T) = -2b 12 (T) / (b n (T) + 
b 22 (T)) , is not strongly affected by temperature for the 
above choice of hadronic parameters: It increases slightly 
with T but remains negative. Hence, changes in the bal- 
ance between the lighter and the heavier species play a 
minor role for the linear approximation. 

Particle densities for nucleons and pions in units of 
n = 0.16 fm~ 3 vs. \i\jm\ = ^ n /m n are shown in 
Figs. 3(a) and (b) for T = 185 MeV. The linear and 
non-linear results are shown for constant excluded vol- 
umes with short dashes and solid lines, respectively, and 
further for relativistic excluded volumes with dotted lines 
and long dashes, respectively. At this high tempera- 
ture the relativistic results are significantly higher than 
the non- relativistic result. A difference between the lin- 
ear and the non-linear approximation due to the non- 
linear enhancement becomes noticeable only for high 
Mn/w n > 0.8 . Thus, for R n = Ri = 0.6 fm from 
above, the linear and non-linear approximation are prac- 
tically equivalent for nucleon densities below n n ss 0.8 no , 

1. e. for densities below about ri\ k> 1/(2 fen). On the 
other hand, due to the strong decrease of b 22 (T) with 
increasing temperature, the influence of the relativistic 
excluded volumes is essential for temperatures of the or- 
der of T ~ TOtt . 

The presence of the additional terms containing tem- 
perature derivatives in the energy density (52) or (53) 
makes it impossible to convert a VdW gas with relativis- 
tic excluded volumes into a gas of free streaming parti- 
cles. Therefore, it is problematic to use these formulae 
for the post-freeze-out stage. For the latter the quantities 
of the free streaming particles without any interaction 
should be used, see discussion in [17,18] and references 
therein. However, these equations of state may be used to 
describe the stage between chemical and thermal freeze- 
out, i. e. a pre-freeze-out stage in terms of Refs. [17,18] . 
This is examplificd in the next section. 
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FIG. 3. Comparison of the model predictions for the nu- 
cleon (a) and pion (b) density n n and , respectively, 
vs. Lin/m n = 0.6 fm, R 2 = 0.3 fm and T = 185 MeV , 
densities in units of no = 0.16 fm -3 ). In both figures the two 
upper lines correspond to relativistic excluded volumes b pq (T) 
and the two lower lines to non-relativistic excluded volumes 
b pq . The results of the linear and non-linear approximation 
coincide - only for extremely large /x n /?rin the non-linear re- 
sults lie slightly higher than the corresponding linear results. 
The deviations due to relativistic excluded volumes are sig- 
nificant. 



V. HARD-CORE RADII 
FROM PARTICLE YIELD RATIOS 

As a simple application of the equations of state pre- 
sented above, let us re-evaluate the thermal model fit pa- 
rameters for particle yield ratios of Ref. [6] , namely the 
hard-core radii of pions and other hadrons R . A two- 
component VdW excluded volume model has been used 
there to explain the pion abundance in A+A-collisions 
by a smaller hard-core radius for the pions than for the 
other hadrons. The ratios has been fitted to BNL AGS 
(Au+Au at 11 AGeV) and CERN SPS (Pb+Pb at 160 
AGeV) data [7] within a thermal model, including all 
resonances up to 2 GeV and using quantum statistics. 

The applied model, however, corresponds to the in- 
correct separated model as pointed out in Sect. I . For 
convenience we give these formulae in Boltzmann approx- 
imation. Within the previously defined notation the two 
coupled transcendental equations read 



?2 P ( T ^l,M2) 



A 1 exp[-(C+£ 2 sp )&ii] 
A 2 exp[-(C+£ 2 sp )M 



(55) 
(56) 



wheras p sp (T, M2) = T (^ sp + £ 2 sp ) • In this context 
A\ represents a sum over the contributions of all hadron 
species but pions, while A 2 corresponds to the pions only. 

The expressions for the particle densities are obtained 
from n| p = dp sp /dp q , 



n 1 " p (T,/ii,^ 2 ) 



n 2 {T,Hi,H2) 



1 + £ P &11+£ P &22 
1+^11+^22 



(57) 
(58) 



Solving these equations for ^ 1 sp and £ 2 sp one recovers the 
canonical pressure formula of the separated model (11) 
as announced in Sect. II . 

Due to the separation of both components in this 
model there is no excluded volume term 612 for the in- 
teraction between different components at all. This is an 
essential difference to both the linear and the non-linear 
approximation. Note that the separated model is not a 
two-component VdW approximation because it cannot 
be obtained by approximating the CPF (13) . 

The transcendental factors of the formulae (55, 56) ex- 
hibit a constant VdW-like suppression exp[— (p/T) b qq ] . 
There is a reduction of this suppression in the linear and 
in the non-linear approximation, as discussed in Sect. Ill . 
The VdW-like suppression is reduced, if 612 appears in 
the corresponding formulae since 612 is smaller than bn 
for R 2 < Ri ■ It is evident that the deviation of the linear 
and non-linear approximation from the separated model 
is the larger the more i?i and R 2 differ from each other. 

In the first step of the fit procedure of Ref. [6] only 
the hadron ratios excluding pions have been taken to 
find the freeze-out parameters. For AGS T w 140 MeV , 
Mb « 590 MeV and for SPS T w 185 MeV , Mb « 270 
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MeV have been obtained. In the second step, a param- 
eter introduced as the pion effective chemical potential 
/i* has been fitted to the pion-to-hadron ratios. Using 
Boltzmann statistics it can be shown that the pion en- 
hancement is thoroughly regulated by the value of fi * [6] ; 
one has obtained [i* w 100 MeV for AGS and fi* w 180 
MeV for SPS data, respectively. 

The pion effective chemical potential depends explic- 
itly on the excluded volumes but also on the pressure. 
The pressure itself is a transcendental function depend- 
ing solely on the excluded volumes since T and jU B 
are already fixed by step one. In Ref. [6] the formula 
H* = (v Q — fjr) p(v , v n ) has been obtained for the sep- 
arated model, where v n = 622 and v = fen are the ex- 
cluded volumes corresponding to the hard-core radii of 
pions R n = R 2 and other hadrons R = R\ , respec- 
tively. Thus, the \i* values for AGS and SPS data define 
two curves in the i? w -i? -plane. The main conclusion of 
Ref. [6] is that the intersection point of these two curves 
(R n = 0.62 fm , R Q = 0.8 fm) gives the correct pair 
of hard-core radii for pions and for the other hadrons, 
i. e. AGS and SPS data are fitted simultaneously within 
the applied model. 

In Ref. [8] these values of R Q and R^ have been crit- 
icised for being unreasonably large. There, a complete 
fit of solely SPS data has been performed within a sep- 
arated model. The best fit has been obtained for equal 
hard-core radii, R v = R Q = 0.3 fm , motivated by nu- 
cleon scattering data. Good agreement has been found 
as well for a baryon hard-core radius, i?Bar = 0.3 fm , and 
a common hard-core radius for all mesons, Rmcs = 0.25 
fm , chooscn in accord with the above ratio of radii, 
Ro/R-n = 0.8/0.62. Larger hard-core radii, especially 
those of Ref. [6] , are quoted as giving distinctly worse 
agreement. 

Assuming the validity of Boltzmann statistics, we have 
re-calculated the i? Q (i? 7r )-curvcs for the above fj,* val- 
ues; firstly in the separated model (55-58) , i. e. as pre- 
sented in [6] . The resulting curves, shown as thin lines in 
Fig. 4 (a) , naturally match the results of the underlying 
fit of Ref. [6] , which arc indicated by markers. 

Then we have considered the linear and the non-linear 
approximation. Due to the occurence of 612 terms in 
these two cases, both functional forms of \i * are essen- 
tially different from the separated case. Consequently, 
the shapes of the i? (i? 7r )-curves are different as well. 
We find distinct deviations from the separated model, 
especially for R^ — ► , and the values for the intersec- 
tion point are slightly lower; see thin lines in Fig. 4 (b) 
for the linear exrapolation. The non-linear approxima- 
tion gives identical results for this purpose because the 
hadron densities arc too small for a noticeable non-linear 
enhancement. 

The crucial point is now to turn on the relativistic 
temperature dependence of the excluded volumes. To 
keep the analysis simple we treat only pions this way 
since they give the strongest effect. 
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FIG. 4. Fits of particle yield ratios for AGS and SPS data 
[7] with the separated model and the linear approximation. 
The thin lines show the fits for the separated model (a) and 
the linear approximation (b) ; the non-linear approximation 
gives identical results for the latter case. 

The thick lines in (a) and (b) show the corresponding curves 
for relativistic excluded volumes v-^ = &22(T) : there is no 
intersection for either of the three models. In both figures 
the results of the fit from Ref. [6] for AGS and SPS data are 
indicated by circles and crosses, respectively. 
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Although the other hadrons are assumed to have equal 
hard-core radii, their relativistic excluded volumes would 
be different for T > , according to their different 
masses. To check the influence of relativistic excluded 
volumes for all particles, we have used one average mass 
of 1 GcV for all other hadrons. The corresponding 
change in the i? -values are below 5% . 

The results of the fit for relativistic excluded volumes 
for pions are shown in Figs. 4 (a) and (b) as thick lines. 
Though this approach is more realistic, there is no inter- 
section point for any of the three models even for very 
large radii R Q , R n ^> 0.5 fm . For the approximated case 
of a single averaged hadron mass there is no intersection 
cither. Because of the different freeze-out temperatures 
for AGS and SPS the v v (T) = 622 (T) values are changed 
differently in both cases, and so are the scales for the 
corresponding R n . 

Due to the errors in experimental data one ought to 
obtain a corridor instead of a curve for each set of data. 
Consequently, the particle yield ratios can be reproduced 
well by e.g. R a w 0.4 fm, R v w 0.2 fm or larger values 
for any of the models with relativistic excluded volumes. 
Therefore, we conclude that the fit procedure proposed in 
Ref. [6] is not suitable to find a unique pair of hard-core 
radii for pions and other hadrons, as long as a best fit is 
searched for just two sets of data of particle yield ratios. 
The use of a relativistic excluded volume for pions along 
with a correct approximation reduce the value of the nec- 
essary nucleon hard-core radius essentially towards more 
realistic values. 



VI. SUMMARY 

In the present work several equations of state for the 
two-component Van der Waals excluded volume model 
are derived and investigated. We have discussed two es- 
sentially different formulations, the linear and the non- 
linear approximation. 

The non-linear approximation is the simplest possibil- 
ity. Here, the large component can reach higher densities 
m than the usual limiting VdW density l/6n, if the 
other component has a suffiently small hard-core radius, 
i?2 < R2, crit • In the linear approximation the densities 
cannot exceed the usual limiting VdW densities I/&11 
and I/622, but generalised excluded volume terms have 
to be introduced. 

For both approximations the suppression factors of the 
grand canonical formulae contain a VdW-like term, pro- 
portional to exp[— (p/T) b qq ] , which however is reduced 
non-trivially. In the linear case there is a slight reduc- 
tion, whcras in the non-linear case this reduction can 
turn the suppression even into an enhancement for the 
smaller component, which leads to exceeding of l/&n for 
the density of the larger component m . 

The commonly used formulae of the separated model 
are shown to be not suitable for the two-component case, 



because they correspond to a system where both compo- 
nents are separated from each other and cannot mix. In 
this model the grand canonical suppression factor is just 
VdW-like and has no reduction of the suppression. 

Furthermore relativistic, i. e. Lorente-contracted, ex- 
cluded volumes have been introduced. Naturally, the rel- 
ativistic excluded volume per particle decreases with ris- 
ing temperature. This effect is the stronger the lighter 
the particle species is. The suppression of particle densi- 
ties in VdW models is lower for a component of smaller 
excluded volume in comparison with a component of 
larger excluded volume. Therefore, the temperature de- 
pendence of the relativistic excluded volumes causes a 
reduction of the particle densities suppression. 

The full equations of state have been presented, for 
both the linear and non-linear approximation, with con- 
stant and with relativistic excluded volumes. For the 
entropy density and energy density there are additional 
terms containing temperature derivatives of the relativis- 
tic excluded volume terms due to their 'thermal com- 
pressibility'. In comparison with the non-relativistic 
case, the expressions for the pressure and the particle 
densities remain unchanged, but the possible range of 
values is obviously wider, since it is generally l/bn(T) > 
I/611 and 1/622(T) > V622 • 

As an application of the derived formulations a fit 
of particle yield ratios for SPS and AGS has been re- 
evaluated. In Ref. [6] this fit had been done in a separated 
model by adjusting the hard-core radii for the pions R n 
and for the other hadrons R . The results of the new fit 
are essentially different from the separated model but co- 
incide for both the linear and non-linear approximation. 
The picture changes drastically, however, if relativistic 
excluded volumes are adopted for pions. The basic idea 
of the fit - one pair of hard-core radii suffices to fit AGS 
and SPS data simultaneously - does not lead to a result 
anymore. This is the case for the separated model and for 
both approximations, linear and non-linear. Experimen- 
tal uncertainties lead to a region of possible values in the 
i?o~-R7i-plane; one could describe the data for R Q > 0.4 
fm and R^ > 0.2 fm . 

We conclude that there are two causes of an enhance- 
ment of particle densities, e. g. the thermal pion abun- 
dance, in VdW descriptions: First, the density suppres- 
sion is generally lower for the smaller component in two- 
component models. Second, there is a further reduc- 
tion of the density suppression due to the relativistic 
excluded volumes. The latter are essentially smaller for 
light hadron species than for heavy species, especially for 
temperatures T > 50 McV . 

When applied to the hadron gas, the linear and non- 
linear results almost coincide for nucleon densities up to 
no ~ 0.16 fm -3 (for R a < 0.6 fm) since the non- linear en- 
hancement does not appear there, but the deviation from 
the incorrect separated model is distinct. However, the 
formulae of the non-linear approximation are essentially 
simpler than these of the linear approximation. 

The influence of relativistic effects on the excluded vol- 
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umes becomes indispensable for temperatures typical for 
heavy ion collisions. Therefore, it is necessary to include 
a correct two- or multi-component VdW approximation 
linear or non-linear - as well as relativistic excluded 
volumes in future calculations. 
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APPENDIX A: THE TWO VDW 
APPROXIMATIONS IN THE CE 

In the following we will study the differences between 
the linear and the non-linear approximation: the total 
excluded volumes v q = v q (Ni, N2) of the corresponding 
partial pressures. In the linear pressure formula (20) each 
component has its own total excluded volume given by 



< nj" < O2) ■ 



(A6) 



,lin 



+ ^2621 , « 2 lin = Nibu + N 2 b 22 , (Al) 



whereas in the non- linear pressure formula (17) there is 
a common total excluded volume for both components 



v 2 



.ni 



^1611+^2622 



JVi N 2 
N ± +N 2 



D . (A2) 



It can be readily checked that it is either t>-/ m < v nl < v 2 n 
or Vi > v nl > v 2 n , i. e. the pole of the non-linear pres- 
sure always lies between both poles of the linear pres- 
sure. Hence, there are values Ni,N 2 , where the non- 
linear pressure is still finite, but the linear pressure for- 
mula is yet invalid since one of the partial pressures has 
already become infinite. Consequently, the domain of the 
non-linear approximation is larger. 

For given V the two domains can be expressed by the 
limiting densities (22) and (23) , which are defined by the 
poles v q (Ni,N 2 ) = V of the corresponding pressure. In 
the linear approximation one obtains the expressions 



1 - 6 2 in 2 „ an, , 1 - 6 22 n 2 

7 , «1,2(«2) = 

Oil 612 

li] 



• (A3) 



For given n 2 , therefore, the domain of p (20) is 

< n{ in < min {n£(n 2 ) , n£(n 2 )} . (A4) 

In the non-linear approximation there is solely one lim- 
iting density 



~ nl/ \_ 1-2 tii 2 n2+\/ (l-26i 2 n2) 2 +46iin 2 (l-fc22"2) / » r\ 
n l \ n 2) — 2fa7I - ^ ' 

Consequently, for given n 2 the domain of p nl (17) is 



In the non-linear approximation there is furthermore a 
region where the pressure has negative partial deriva- 
tives with respect to the smaller particles' number, 
dp al /dN 2 < 0. The condition dp al /dN 2 = defines 
the boundary of this region 



(A7) 



- nl. bd / \ 

n x (n 2 ) 

__ 1-2 (612-622) n 2 +- v /(l-2 (b 12 -b 22 ) k 2 ) 2 +8 (611-612) n 2 
4 (611-612) 

For given n 2 a negative derivative dp nl /dN 2 < oc- 
curs only at a density n\ > n " ' bd , while the deriva- 
tive dp nl /dNi is always positive for R 2 < Ri as readily 
checked. 

In Fig. 1 (a) the functions rii(n 2 ) are presented in di- 
mensionless variables ni&n and 712622 ■ The properties 
of these dimensionless functions depend only on the ra- 
tio of the two radii R 2 /R\. The smaller this ratio is, 
the higher is the maximum value of n" 1 , while the re- 
gion of negative derivatives dp nl /dN 2 becomes narrower. 
The straight line h^\(n 2 ) starts always at I/&11, but 
its slope decreases for smaller R 2 /R\ , whereas h 2 {n 2 ) 
ends at I/622 and its slope increases. The pressure of the 
separated model (11) would yield one straight line from 
ni&u = 1 to n 2 6 2 2 — 1 in Fig. 1 (a) , for any ratio of the 
radii i?i and R 2 . 

For very small ratios R 2 /R\ , i. e. for R 2 — > , one finds 
from Eq. (A2) that v nl -» Nrfu [1 - \N 2 /{N X + N 2 )} . 
This yields the maximum density max(n 1 nl ) = 4/6n for 
N 2 3> Ni . Thus nf 1 exceeds the maximum density of 
the linear approximation or of the corresponding one- 
component VdW gas, max(nj™) = max(n 1 )C ) = l/6n , 
by a factor of four in this case. 

Note that the value 4/&n appears in the linear approx- 
imation as well: For 7j 2 lln — > V it is max(nj") = 4/6n at 
n 2 = , but this density cannot be achieved because 
Pi n (ni,n 2 ) is infinite for m > l/6n . 

Let us consider now the consequences of negative 
derivatives dp nl /dN 2 < in the non- linear approxi- 
mation. If a negative dp nl /dN 2 occurs for a density 
n{ = const, at n 2 = 0, the pressure p nl (n[,n 2 ) has 
a minimum at a certain density n 2i m i n > , which is 
determined by the boundary n^ l ' hd (n 2 ) . For increas- 
ing ni along the boundary, consequently, the non-linear 
pressure is always lower than for increasing m at fixed 
n 2 = . Hence along the boundary higher densities can 
be achieved, in particular ni > l/6n . 

Therefore, exceeding of n" 1 = l/6n requires that the 
boundary starts inside the the non-linear domain at n 2 = 



. Thus the condition 
critical radius R 2 , cr i t (24) , 



- nl,bd 



(0) < l/6n provides the 



611 < 2612 



#2,crit(i?l) = (v^-l)iZi . (A8) 



On the other hand the boundary starts at ni = 8/ (14 bu) 
for R 2 — > at n 2 = , i. c. for any density 6nni < 



13 



8/14 w 1/2 negative values of dp nl /dN 2 do not occur for 
any radii. 

Although it is pathological that for high densities n\ 
the non- linear pressure firstly decreases, if particles of the 
second and smaller component are added to the system, 
there is a reasonable explanation for the lowered pressure 
along the boundary (A7) for small radii R2 < R2, crit • 

Consider, for instance, n\bu — 0.9 in Fig. 1 (a) . Since 
it is R2/R1 — 0.4 , the dimensionless density of the small 
particles at the boundary n" 1,bd (n2) nearly vanishes, 
ri2&22 ~ 0-05 7 whereas the absolute amounts of the small 
and large particles are about equal. For the excluded 
volume interaction of the large particles in the pressure 
formula (17) , therefore, the influence of the mixed term 
612 becomes comparable to that of the distinctly larger 
non-mixed term bu . 

For R2/R1 <C 1 one obtains n 2 ^> n\ near the bound- 
ary at m&n = 0.9 , i. e. the large particles are completely 
surrounded by the smaller particles and interact mostly 
with these but hardly with other large particles anymore. 
In this situation, consequently, the excluded volume in- 
teraction of the large particles is governed by the essen- 
tially smaller 612 and not by bu < 8 612 . 

One might interprete this behaviour as an effective at- 
traction between small and large particles, but it is rather 
a strong reduction of the large particles' excluded volume 
suppression. 

As VdW approximations are low density approxima- 
tions they coincide for these densities, but they evidently 
become inadequate near the limiting densities: both dis- 
cussed formulations do evidently not match the real gas 
of rigid spheres there. 

For high densities the linear approximation behaves 
natural, i.e. it is dp hn /dN q > always. However, one 
has to introduce the additional terms 612 and 621 . For 
the choice (21) these terms provide a one-component-like 
behaviour in the limits R2 = Ri and R2 = , but they 
have no concrete physical meaning. 

In the non-linear approximation there occur pathologic 
pressure derivatives dp nl /dN2 < for R2 -C i?i . How- 
ever, the non-linear formulae may be used for special pur- 
poses, e.g. for m > I/611 at intermediate n 2 &22 , where 
the linear approximation is yet invalid. 



Mp (T, V, N U N 2 ) = —T qjj- ln[Z(T, V, N,, N 2 )} . (Bl) 

Partial derivatives of the pressure with respect to the 
particle numbers dp/ dN q have no relevance here. 

For the examination of chemical stability it is ap- 
propriate to switch from the free energy of the CE, 
F(T,V,N 1 ,N 2 ) = -T ]n[Z(T,V,N 1 ,N 2 )] , to the Gibbs 
free energy or free enthalpie 



G(T,p,N u N 2 ) = F+pV = n 1 N 1 +ii 2 N 2 



(B2) 



where fj, q (T,p, N 1: N 2 ) = dG/dN q . This requires that 
p(T, V, N x , N 2 ) can be solved for V(T,p, N U N 2 ) , which 
is the case for the non-linear approximation, 

V n \T,p, N l7 N 2 ) = + iVi&n + N 2 b 22 - w^k D ■ 

Further it is useful to introduce the molar free enthalpie 
q = G/(Ni + N 2 ) — g(T,p,x\) with the molar fractions 
xi = JVi/(JVi + N 2 ) and (1 - an) = x 2 = N 2 /(N 1 + N 2 ) 
of component 1 and 2 , respectively. Then the chemi- 
cal stability of a binary mixture [19] corresponds to the 
condition 



dx- 



■g{T,p, Xl ) 



dx± 



dii 2 {T,p,x 1 ) 
dxi 



> 



(B3) 



For the non-linear approximation one obtains 
3 nl (T,p,xi) 



Xl 

+ 



(l"*l){ 



P (bu 



T In 



l-xi P 
<t>2 T 



(1-Zl) 

+ P (b 2 2 



D)} (B4) 
■x*D)} , 



and thus condition (B3) is satisfied: 

^Lg»\T,p, Xl ) = i + ^ + P 2D > . (B5) 



Therefore, the system described by the non-linear ap- 
proximation is thermodynamically stable - despite the 
pathologic behaviour in special cases. Due to the equiv- 
alence of the thermodynamical ensembles this is true for 
any representation of the model. 



APPENDIX B: STABILITY OF THE 
NON-LINEAR APPROXIMATION 

The non-linear enhancement in the GCE or the oc- 
curence of negative values for dp nl /dN q in the CE sug- 
gest a further investigation concerning the thermody- 
namical stability of the non-linear approximation. 

One can readily check that in the CE it is dp nl /dV < 
generally, so there is no mechanical instability. To in- 
vestigate wether there is a chemical instability [19] it is 
necessary to study partial derivatives with respect to the 
particle numbers, d/dN q , of the chemical potentials 



APPENDIX C: THE TWO VDW 
APPROXIMATIONS IN THE GCE 

In this part we will study the non-linear and the linear 
approximation in the grand canonical ensemble. As in 
the CE the differences between the linear and non-linear 
approximation occur only for high densities of the larger 
particles, n q , which correspond to large chemical poten- 
tials [i q in the grand canonical treatment. Therefore, we 
will study the limit given by Eq. (44) : \i\jT — ► 00 for 
constant T, /12 and R2 < Ri , where it is £1 — > 00 . 

The transcendental exponents of both approximations 
contain an attractive part besides the usual VdW-like 
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suppressive part exp[— (p/T) : In the linear approxi- 
mation (36, 37) the suppression is reduced, but the tran- 
scendental exponents are always negative - whereas the 
suppression in the non-linear approximation (29, 30) is 
not only reduced, but the transcendental exponent of £ 2 nl 
can even become positive in the above limit. 

To examine the latter we rewrite the coupled transcen- 
dental equations (29, 30) as 



d nl = 0i exp [£] x 



x exp 



- (?i nl + e 2 nl ) &n i- 



D/bi 



£ 2 nl = 2 exp [f ] x 



(CI) 



(C2) 



x exp 



(Ci nl +e 2 nl ) & 22 1 



D/b 2 



(i+^V?" 1 )' 



If R 2 is sufficiently smaller than R\ , then D/b 22 is larger 
than unity and the transcendental exponent of £ 2 nl (C2) 
can become positive, 



< - (\ + ^ 



«i nl 



< a 2 



(C3) 



The coefficient a 2 = a 2 (Ri/R 2 ) introduced here is the 
crucial combination of excluded volumes in the non-linear 
approximation. It characterises the behaviour of this ap- 
proximation in the limit (44) , i. c. for £1 — > 00 . 

For equal radii R 2 = Ri it is a 2 = — 1 , and one has full 
VdW-like suppression, cx cxp[— (p/T) b qq ] . Negative a 2 
indicate the strength of suppression of £ 2 nl for increasing 
. For — 1 < a 2 < the suppression is reduced, most 
strongly for a 2 w . In the case a 2 — there is no 
suppression in the limit (44) but £ 2 nl —fA 2 = const, for 
— ► 00 . Thus, the condition a 2 = provides the 
critical radius, 



D/b 22 = 1 



R 



2, crit 



= (&4-l)R 1 



(C4) 



which coincides with the canonical result (A8) . 

For positive a 2 or R 2 < #2, crit the non-linear en- 
hancement occurs: £ 2 nl is enhanced by increasing as 
long as the second exponent in (C2) is positive, i. e. for 
£ 2 nl < a 2 d nl . The transcendental factor of (CI) has 
only a reduced suppression in this case. According to 
Eq. (C3) one obtains 



i. 



ni 



«2 i 



ni 



but also 



a 2 n x 



ni 



(C5) 



since it is nfjnf = £ 2 nl /£i nl due to Eqs. (33, 34) . 

Using Eqs. (C5) and (C3) one obtains for the particle 
densities (33, 34) 



nf 



bn-ag b 22 fcn-Cv^D-V^) 2 



a 2 n x 



ni 



a 2 _ 



Vd- 



bu-a* b 22 622 6n-(V^-\/E22 ) 2 



(C6) 



(C7) 



It is clearly seen that the density n" 1 can exceed I/&11 
for positive a 2 . As in the CE the maximum value, 
max(n"') = 4/6n , is achieved for R 2 = or a 2 = 00 . 
Then, in the limit (44) the density of the second compo- 
nent diverges, n 2 x — > 00 , but its total excluded volume 
vanishes, n 2 x b 22 — > , as seen from Eq. (C7) . 

There is yet another case, where the condition (C3) 
is fulfilled, the early enhancement : £ 2 nl can be enhanced 
with increasing /i 2 for constant T and n\ , if fii is suffi- 
ciently large. This enhancement takes place only at small 
[i 2 , and it obviously vanishes when \i 2 is large enough so 
that £ 2 nl > a 2 £ 1 nl . The early enhancement is the direct 
analogue to a negative derivative dp nl /dN q < in the 
CE. 

The coupled transcendental equations of the linear 
approximation (36, 37) may be rewritten similarly to 
Eqs. (CI, C2) . For the choice (21) one obtains in terms 
of D from Eq. (16) 



ercr,Mi,M2) 



= A\ exp 



£ 2 lm (T, m , M2 ) 
= ^4 2 exp 



(« 



(C8) 

)]■ 

(C9) 



In this case the condition for the enhancement of £ 2 
R 2 < Ri would be 



lin for 



< a 2 



~ _ 2612 



611+622 



(CIO) 



As a 2 is always negative the VdW-like suppression is 
only reduced in this approximation. Like in the non- 
linear approximation a 2 — — 1 corresponds to equal radii 
R 2 = Ri and full VdW-like suppresion, whereas a 2 = 
— 1/4 (R 2 = 0) corresponds to the most strongly reduced 
suppresion in the linear approximation. 

Thus, the densities n q ln can not exceed the maximum 
value 1/bgq of the corresponding one-component case. 
Furthermore, the density n 2 ln (41) vanishes in the anal- 



d 1 



lin 



00 , in contrast to the non- 



ogous limit to (44) . 
linear behaviour. 

In the case R\ < R 2 one w ould evidently investigate 
the coefficients a\ = (y/D/bn — 1) in the non-linear and 
ai = — (26i 2 )/(6n + b 22 ) = a 2 in the linear approxima- 
tion, respectively. 
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